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We establish the relation between the Wigner- Weisskopf theory for the description of an unstable 
system and the theory of coupling to an environment. According to the Wigner- Weisskopf general 
approach, even within the pole approximation the evolution of a total system subspace is not 
an exact semigroup for multichannel decay, unless the projectors into eigenstates of the reduced 
evolution generator W{z) are orthogonal. With multichannel decay, the projectors must be evaluated 
at different pole locations Za 7^ zp, and since the orthogonality relation does not generally hold 
at different values of z, the semigroup evolution is a poor approximation for the multi-channel 
decay, even for very weak coupling. Nevertheless, if the theory is generalized to take into account 
interactions with an environment, one can ensure orthogonality of the W{z) projectors regardless 
the number of the poles. Such a possibility occurs when W{z), and hence its eigenvectors, are 
independent of z, which corresponds to the Markovian limit of the coupling to the continuum 
spectrum. 
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PACS numbers: 

INTRODUCTION 

Many physical systems demonstrate instability, i.e., a 
transition from a relatively stable state to a final state 
which in general corresponds to a system with many iden- 
tifiable degrees of freedom. Such occurs, for example, in 
particle decay or radiative atomic transitions. In many 
cases the process observed has the semigroup property, 
i.e., the operators generating the evolution on the Hilbert 
space of quantum states satisfy the composition law 



Z{ti)Z{t2)^Z{ti+t2) 



(1) 



for ti,<2 > 0, and for which the operators Z{t) do not 
have an inverse and are decreasing. This relation is a 
general form of the well-known exponential decay law, 
for which Gamow [l[ constructed a phenomenological 
Schrodinger equation with complex energy eigenvalue of 
negative imaginary part. Weisskopf and Wigner Q for- 
mulated a basic theory which approximately reproduced 
the Gamow result in second order perturbation for a sin- 
gle channel decay. In their original formulation, the sur- 
vival amplitude of a quantum state \ip) is given as 



-iHt 



(2) 



where H is the full Hamiltonian of a system consisting of 
an unperturbed part Hq for which ji/)) is an eigenstatc, 
and a perturbation V is understood to induce a transi- 
tion to an infinite number of final states with a continu- 
ous spectrum. The Laplace transform of Eq. ([2|) provides 
an expression corresponding to the Green's function (or 
resolvent kernel) for the Schrodinger evolution. In a sin- 
gle channel decay problem the pole approximation for 



the inverse Laplace transform results in an approximate 
semigroup property for U^'^'^{t), which is equivalent to 
the perturbative analysis of Weisskopf and Wigner. How- 
ever, in case of the two or more channel decay, such as the 
neutral K meson decay (for which there is CP symnmetry 
breaking) , it has been shown that the pole approximation 
does not reproduce the semigroup evolution observed in 
Q. To account for semigroup behavior application has 
been made of a quantum forniH, @| of the classical Lax- 
Phillips theory 16[ (see also [7[ and references therein) 
and its generalizations @ . This theory achieves an exact 
semigroup law by imbedding the usual quantum theory in 
a larger Hilbert space, consisting of a direct integral of a 
family of Hilbert spaces of usual type, foliated according 
to the time parameter. The result was particularly effec- 
tive and straightforward for treating quantum mechan- 
ical systems with Hamiltonians of unbounded spectrum 
llCf . Its generalization for problems with semibounded 
spectrum Q made it clear that the imbedding associated 
with the Lax-Phillips theory effectively introduces many 
additional degrees of freedom. In this paper, we show 
that semigroup evolution may be obtained working in 
the framework of the usual quantum theory through an 
explicit coupling to many environmental degrees of free- 
dom. Thus, we provide a physical framework accounting 
for the mathematical structure of the Lax-Phillips theory. 

We study the Wigner- Weisskopf pole approximation 
theory Q in the framework of the Lcc-Fricdrichs model 
This model consists of a subspace of discrete 



states interacting with a continuum subspace of states, 
for which there is no direct continuum-continuum 
interaction. It is convenient to take the discrete-discrete 
interaction to vanish as well. The association of the 
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continuum with an environment constitutes a new aspect 
of the model. The "environment" here may be thought 
of as a distribution of a very large, or even infinite 
number of final states into which the initial state of the 
system decays. The Lee-Friedrichs model ll| in Lee's 



construction was formulated in terms of non-relativistic 
quantum field theory using a Hamiltonian for which 
the interaction l/^ee ^ affayOQ + h.c. included an 
annihilation operator ajv for the original (unstable) state 
N, multiplied by creation operators ayOg for the two 
body final state V and Q. These operators could have 
been constructed to include creation of a many-body 
environment as well, maintaining its equivalence with 
the Friedrichs quantum mechanical form, with spectral 
coefficients coupling to a bath. In other words, the 
decay process may include not just the specified final 
states of the decay model, but also an environment. It is 
instructive to think of the spontaneous emission process 
12, [13] as the most well-known illustration of a decay 
into an infinity of final states, where an excited atomic 
state decays into a distribution of Fock space states of 
the radiation field (i.e. photons). 

In what follows we show that by associating the 
continuous spectrum with an environment the reduced 
evolution of the discrete states subspace is defined by 
a spectral correlation matrix a{t), which is a function 
unifying the details of the interaction. The notion 
of the spectral correlation matrix is inspired by the 
well-studied particle(8'environment theories [13] ■ In these 
models the reduced evolution of the particle is obtained 
by tracing out the environmental degrees of freedom, 
leaving a stochastic dynamical equation which invokes 
the so-called environmental (in some examples complex) 
noise z{t). As shown in [15|, the environmental correla- 
tion function corresponds to the autocorrelation value 
of this noise, i.e. a{t) = {z*{t)z{0)), and in agreement 
with the fluctuation-dissipation theorem turns out to be 
a time dependent memory-kernel for the particle energy 
dissipation. Exploiting the analogous properties of the 
environmental correlation function and the spectral 
correlation matrix we investigate the validity of the 
non-trivial semigroup evolution for the many channel 
decay and identify the necessary conditions with the 
well-known Markovian limit [l2|, [ij]. In this way we 



consider a natural imbedding of the Wigner-Weisskopf 
idea into a theory of interaction with a reservoir (e.g., 
generalization of Anderson and Fano 16| and Lee's 
formulation) . 

The next two sections include a review of essential 
results from the Wigner-Weisskopf theory and the 
derivation of the Markovian limit, provided with a brief 
summary of situations where this limit may be realized 
exactly or approximately (readers familiar with these 
concepts are welcome to quickly leaf through). In the 



fourth section, which includes the main novelty of the 
paper, we discuss the effect of the Markovian limit on 
the Wigner-Weisskopf pole approximation, and explain 
how the semigroup evolution law is achieved for the 
many channel decay as well. 



THE WIGNER-WEISSKOPF METHOD 

In this section we summarize a functional formulation 
of the Wigner-Weisskopf pole approximation theory [^ 
and its generalization to many channel decay [3] . Follow- 
ing the usual model for the decay of an unstable system 
we consider a Hamiltonian of the form 



H = Ho + V, 



(3) 



where the spectrum of Hq consists of a finite number 
N of discrete eigenvalues {Aq}, embedded into a con- 
tinuum {A > 0} with spectral weight dE{\) = \\)(X\dX. 
We study, in particular, the Lee-Friedrichs model [ill ], 
for which the interaction V couples the discrete states 
to the continuum, but does not couple continuum states 
or discrete states among themselves. The fact that the 
continuum subspace is associated with the products of 
the decay process presents an opportunity for the intro- 
duction of an environment. We assume that the initial 
unstable state of the system is given by a superposition 
of the discrete eigenstates {]</>«)} oi Ho: 



N 

E 



Ca\(l}a); Ho\(t)a) = \a\4>a)- 



(4) 



The reduced evolution, i.e. the evolution of the dis- 
crete subspace, is governed by the reduced propagator 
R{z). This propagator in Laplace transform is defined 
(for Im z > 0) by projection of the total system propa- 
gator 



U{z)= / dte'^'e 



t„-iHt 



z-H 



(5) 



on the discrete subspace. The a/3 matrix element of the 
reduced propagator R{z) in the Laplace domain is 



Rap{z) = -i((f>a\U{z)\(l)p) = {(f>a\ 



1 



Z-H' 



3)- (6) 



Eq. ^ may be written using N x N matrix notation 
(defining W{z)) 



R{z) 



z~w{zy 



(7) 



which is confined to the discrete subspace, and where, 
by comparing with Eq. W{z) is understood as the 
Laplace space reduced evolution generator. The reduced 
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evolution in the time domain, dictated by U'^'^'^{t), 
given by the inverse Laplace transform of Eq. ([7]) 



1 



R{z)e''''*dz, 



c 



where the contour of the integration C, shown in Fig. [1] 
runs slightly above the real line on the z-plane from +oo 
to zero (the bottom of the positive spectrum), and then, 
around the branch point, from zero back to +oo slightly 
below the real line. 



described in detail in Q , results in the following alterna- 
tive expression for the reduced propagator 



(8) rr 



Res 



(10) 

where Res[R"{zj))] {R" {z) on the right and i?(z), 
which we denote above by R^{z), on the left) is the 
residue of R^^{z) at the pole position Zj and the contour 
of the integration Ci , shown in Fig. [21 now runs around 
the branch point along the negative imaginary axis. 
The integration along the contour Ci carries the factor 
g-!zt 2 ^[jg lower half plane. Hence, for i > and 
not too small^ one can consider neglecting this term, 
called the "background contribution". (There is a very 
long time contribution from the neighborhood of the 
branch point which we do not consider here). Doing so, 
the evaluation of the reduced propagator reduces to the 
summation of the contributions of the residues of the 
poles of R^^{z) in the lower half plane. The assumed 
dominance of the these contributions is called the "pole 
approximation" . 



FIG. 1: The inverse Laplace transform contour C in Eq. 



Since, in the general case, the exact calculation of 
Eq. ([8]) is difhcult, one is interested in a useful approxima- 
tion. We first note that there can be no pole for Im z ^ 
in Eq. ([7|) on the first Riemann sheet (see Appendix). 
However, we may explicitly continue the integration in 
Eq. ([5]) analytically to the second Riemann sheet (using 
Eq. (HI)) 0.' Doing so yields 



i?"(z) = 



(9) 



z- VF"(2)' 

which allows one to deform the contour of integration 



To obtain an expression describing the residues of 
R^^{z), we focus on the reduced generator W^^{z). Even 
though W^'(z) is not generally self adjoint, Eq. ^ may 
be represented in a way analogous to the spectral the- 
orem, as the sum of normalized projectors {Qaiz)}, so 
that 



TV 



where 



Qa{z) 



\a,z)R L{a,z\ 
L{a, z\a, z)r 



(11) 



(12) 



are made of the left and right eigenvectors of W^^{z) 
(corresponding to appropriate linear combinations of 
the eigenvectors of Hq, sometimes called "decay eigen- 
states"): 



W"(z)|a, z)r = uja{z)\a, z)r, 



,(a, z\rW^^{z) = u:a{z)L{a, z\. 



(13) 



FIG. 2: The modified inverse Laplace transform contour Ci 
in Eq. (|10[) . The poles are denoted schematically by the stars. 

The eigenvalues of VF^^(z) (see below) determine the 
poles of R^^{z) in the lower half plane. The procedure. 



[1] The scale of the value necessary to go beyond the non-exponential 
region of decay curve, often called the Zeno time [TtI , is deter- 
mined by the dispersion of the Hamiltonian, as discussed in [2lll 
and below (Eqs. (37), (38)). The onset of the exponential behav- 
ior generally occurs after the curve of steepest descent, which 
rotates clockwise in time [2ll |. passes the first pole, which then 
dominates the time dependence. 
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Clearly, 

z|W^"(z)|a, = uja{z)L{l3,z\a,z)R = 

= UJi3{z)L{l3,z\a,z)B. (14) 

may be valid for iOa{z) ^ ^^Jpiz) (for any z) only if 

i(/3,z|a,z)i? = 0. (15) 

This orthogonality relation for the eigenvectors of W^^{z) 
provides the orthogonality of the appropriately normal- 
ized projectors 

Qa{z)Qp{z) = Qa{z)SaP, (16) 

at each point z. Eq. PT|) follows from Eq. © 
by using the spectral representation Vl^^^(z) = 
J2a '-^a{z)\a, z) fi L{a, z\ and the orthogonality proper- 
ties of the Qa{z)- 

Eq. (jl6p plays crucial role in examination of the semi- 
group property of the reduced evolution. Applying the 
pole approximation procedure, we neglect the "back- 
ground contribution" and approximate the reduced prop- 
agator by the sum of the residues of R^^{z) Eq. 
which for weak coupling V may be well approximated 
to yield 



ed 



j 



^aAZj), 



(17) 



where aj corresponds to the singularity at Zj — ujaj{zj). 
Repeated application of this reduced evolution is then 



e-"^''e-'^''''Qo.^iz,)Q^,{zk). 



U'''"^{t2W'^{ti)^Yl 

. 

Although for the single channel problem, if there is just 
one pole, the projectors product in the right hand-side 
of the last equation is trivially unity and Eq. (jl8p shows 
semigroup decay, for the many channel decay with many 
poles the projectors Qaj{zj),Qak{zk) are generally eval- 
uated at different pole locations on the Laplace plane. 
For Zj 7^ Zfc the orthogonality relation Eq. ((T5)) can no 
longer ensure Qaj{zj)Qah{^k) = 0. Thus, even in the 
pole approximation, the semigroup evolution is generally 
not valid ISHlQi, i.e. 



-iZj (tl+t2)/ 



(Zj). 
. 

In spite of this conclusion, many experiments display 
semigroup decay to high accuracy [3| , while the estimates 
lil lij have shown that the deviations predicted by the 
Wigner-Weisskopf theory would exceed the experimental 
error 



ASSOCIATION OF THE SPECTRAL DENSITY 
FUNCTION WITH AN ENVIRONMENT 

In this section we sketch a treatment for the reduced 
system dynamics familiar in the field of quantum op- 
tics 12 , 3] and condensed matter physics [1^ [l6| . We 
adopt the Lee-Friedrichs model of the system described 
in the previous section for the total Hamiltonian given 
by Eq. ([3]). For our current purpose it is not necessary to 
make any preliminary assumptions (other than, for sim- 
plicity, degeneracy) regarding the structure of the contin- 
uous part of Hq] it may be bounded from below or may 
not. The interaction picture propagator, defined through 
the total system propagator U(t2 — ti) and the unper- 
turbed total system propagator Uo{t — to) = Q-iHa(t-to) 
is given by 

U{t2M)=U^\t2~tQ)U{t2~ti)UQ{ti~to) (20) 



and obeys the equation 

. d ~, 

i-rU{tM) 
at 

where 



= Vlt)lHtM), 



Vit) = U^\t^to)VUo{t-to) 



(21) 



(22) 



is the interaction picture Hamiltonian. Integrating 
Eq. pip and iterating it one time we get the exact equa- 
tion 

U{t,to)^l-if V{T)dT-f I V{T)V{T')U{T',to)dT'dT. 
J to J to J to 

(23) 

Projecting the last expression on the discrete subspace 
we find 

U'^f{t,to) = {^a\U{t,toMp)^ 



E 



to J to J \ 



iA„(r-to) 



{<Pc.\V\X)e 



'iX(T-ta) , 



X e^^^-''-^°\\\V\(t)^)e-'^-'^^'-^°'^lJ:^f{T'M)d\dT'dT. 

(24) 

Here the first order term {(f)a\V{T)\(l)fi) has vanished, be- 
cause by assumption, V does not couple the discrete 
states among themselves, and the last term was obtained 
using Eq. ^F^. Next we differentiate Eq. ([^ with re- 
spect to t and obtain the reduced, i.e., projected into the 
unstable subspace, master equation 



-^Uaf'{t,to) 



. y^gzAo(t-to). 
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J to J \ 

(25) 

where the matrix elements 



w,^(A) = {(^^\V\\){\\V\<P-, 



(26) 



form the so-called spectral density matrix. The integral 
transform of the latter 



(27) 

defines the elements of the spectral correlation matrix 
within the interaction representation. Using the last ex- 
pression, Eq. ([25]) may be written as 



dt 



KfitM^^Y. f 5ca.,{ty)U;f{T'M)dT'. (28) 



Transforming back to the Schrodinger representation, we 
use the inverse version of Eq. (j20p . Taking into account 
that the spectral correlation matrix a{t,T') transforms 
analogously to U{t2,ti), that JjC/(t,io) — iHoUQ^{t - 
to)U{t,to) + UQ\t - to)j^U{t,to), and that Hq and 
^^(^2 — ^i) are diagonal in the reduced basis {|0c<)} rep- 
resentation, in terms of matrix notation we obtain 



U'^'^^'it, to) - -iH^'^U'-'^'it, to)- a{t-T')U^^^iT', to)dT'. 

■J to 

(29) 



dt 
Here 



N 



(30) 



is the discrete part of the unperturbed Hamiltonian Ho, 
and the spectral correlation matrix defined as 



ait) 



-i\t 



w(A) d\ 



(31) 



is the Fourier transform of the spectral density matrix 

N 



w(A) = ^ (A)|0a)(<^, 
Q,7— 1 



We argue that a{t) defined by Eq. (|3T|) can be as- 
sociated with the "noisy" environmental correlation 
function ait) = {z*it)ziQ)) mentioned earlier, since the 
microscopic definition IJ, |l5| of the latter evidently 
coincides with Eq. (|3T|) up to an obvious generalization. 
The correlation matrix ait) Eq. ([31]) represents all the 
microscopic details of the interaction, whose properties 
determine the type of the reduced evolution, as will be 
clear from the following. 



Equivalence of Exact Markovian Coupling and 
Globally Flat and Unbounded Spectral Density 
Matrix 

Note, that in an ideal case for which LOa-yiX) Eq. (pS)) 
are independent of A and the continuous spectrum of 
Ho is unbounded (a physical example of such a situ- 
ation could occur for a Stark type interaction with a 
bath, which induces a shot noise), the correlation ma- 
trix Eq. ([5T|) reduces to a delta function of time: 



iX{t- 



^>d\ = 2TTUj5it - t). (33) 



ait 



Substituting this into Eq. (|29|) we find 

jU'^'it - to) = i-iH^"" - 2t:u) U^-'^it - to). (34) 

The last equation is a local in time first order differen- 
tial equation with constant evolution generator, i.e. an 
equation describing semigroup evolution, also called the 
Markov equation, since it describes a Markovian stochas- 
tic process [23|. Its solution for t > is 



(35) 



where the decay matrix T is given by 

r E 



t .t 

aiT)dT = / 2'KUj5it - T')dT = 2ttoj. (36) 

to ^to 



The well-known demonstration of the Zeno effect for very 
short times t relies on an expansion of the survival prob- 
ability Pit) of a unstable state \ip) in a series 211 



Pit) = |(V'|e-'"*|V') 



- iHt - -H^t^ + c^il- t'^AH^ (37) 



where 



AH' 



(38) 



is the dispersion of the total Hamiltonian H in the state 
(32) The standard argument leading one to a conclusion 

about the principle impossibility of pure semigroup evo- 
lution is the presumable possibility to cut off the above 
expansion after the second order in t. Note, however, 
that for the Lce-Fricdrichs model in the Markovian limit 
the coefficient AH^ diverges. If, for example, \tp) — \ipa), 
Eq. dSS]) becomes 



/OO 
dXojU^) ^ 00 
-00 

(39) 

Hence, the truncation of the expansion is invalid together 
with the physically interpreted conclusion. It should 
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be stressed, that under our assumptions regarding the 
structure of the coupling, there is no Zcno (sec also (22| ) 
effect and the result Eq. ([55]) is exact. 



Now wc briefly review how such assumptions may be 
realized. For this purpose we adopt the usual approach to 
the original definition of the spectral density function, i.e. 
as the continuum limit of a quasi-continuous spectrum. 
We assume, that the unperturbed system, described by 
Hq, is confined to a large box with some standard bound- 
ary conditions. The spectrum of Hq then consists of 
A'' discrete eigenstates embedded into a quasi-continuous 
spectrum of the environment including the decay prod- 
ucts. Each element a;Q.^(A) of the spectral density matrix 
a; (A) Eq. (f32|) is defined as the product of the interaction 
amplitudes {(pa\V\X) {X\V\4>^) , weighted by the local de- 
generacy of the environmental states 



UJajW = {(bo.\V\X){X\V\cl)^)D{X), 



(40) 



where the density of states D{X) serves as an effective 
"coarse graining" . When the spectrum is truly continu- 
ous D{X) may be absorbed in {(l)a\V\X) {X\V\(j)^)dX. The 
interaction amplitudes ((/)Q|y|A)(A|y|0^) are determined 
by the microscopic details, so their functional depen- 
dence on A is dictated by the nature of the interaction. 
On the other hand, the density of the environmental 
states D[X) is determined by Hq and, in particular, 
by the boundary conditions to which the system is 
confined. Thus, D{X) is at our disposal. Designing 
the geometry of the space we may find a D(X) which 
suppresses, at least approximately, the A-dcpendcncc 
of ((/>Q,|y|A){A|V'|</i-).), and yields a semigroup evolution 
with a desirable accuracy^. 



Approximate Markovian Coupling vs. Weak 
Interaction and Resonance. 



In the preceding subsection we have assumed a contin- 
uous spectrum on the whole real line in order to obtain 
an exact semigroup evolution. However, if there is a res- 
onance and the coupling is weak, the leading behavior of 



[2] There exists a class of phenomenological descriptions of the spec- 
tral density functions aJa'y(A) oc rysAc"" A''e~^''^'= []3 |. where 
ris is the viscosity constant, the exponential factor provides 
a smooth cut-off modulated at the frequency A^, and wrhcrc 
< s < 1 and s > 1 describe the so-called sub-ohmic and the 
super-ohmic interaction respectively. The slower is the depen- 
dence of the spectral density function on A, the closer is the 
reduced evolution to the Markovian limit. The boundary case of 
s = 1 corresponds to the ohmic interaction, as for example, the 
dipole interaction of a particle with the free electromagnetic radi- 
ation field, which is still associated with an approximate Marko- 
vian spectral coupling. 



the system is dictated by small regions of the continuous 
spectrum where the effect of the interaction is sharply 
enhanced. Such regions occur at A sufficiently close to 
the resonances, which for small coupling are close to the 
eigenvalues {Xa} of the discrete subspace. To show this 
we return to Eq. (^5]) and see whether it is possible to 
utilize it approximately, even if w(A) is not constant and 
the continuous part of the Hq spectrum is bounded from 
below. Inspecting the integrand of the last term on the 
right hand-side of Eq. ([25|) we note that the collapse of 
the memory kernel may result from rather symmetrical 
manipulations with respect to the spectral or the time 
variables. For the perfectly Markovian coupling a; (A) 
is constant and unbounded, the integral over A yields 
S{t — t') and admits the semigroup evolution of Eq. ([M]) . 
Conversely, if it is justified to neglect the difference be- 
tween {/^^''(r', io) and U^'jfit, ta) and to stretch the hmits 
of the integration over t' to ±oo, it also yields a Marko- 
vian equation of the form of Eq. (|33p (and results in a 
factor (5(A — A^) as well). The error caused by the replace- 
ment U:^f{T',to) UJ^fit.to) (X 0{V^) is negligible in 
case the interaction is sufficiently weak, because the last 
term of Eq. (pS)) is already second-order in V . Such an 
approximation yields a solution, which up to the second 
order in V is exact, and provides a direct interpretation of 
the original perturbative computation of Wcisskopf and 
Wigner 0] in the resolvent formalism. Making the weak 
coupling assumption, wc approximate the reduced master 
equation (^5]) by 



E 
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giA„(t-to)giA^tog- 



iXt, 



(A)) 



^i^-,~>^)^'dT'dx tj':f(tM) 



(41) 



and focus on the integration over the time. Since the 
interaction V is time independent, one may select the 
time origin such that to = —T/2 and t = T/2. Then the 
integration over r' in Eq. (|^T|) yields [l2| 



T/2 



T/2 



-^^^'dr' - - 

TT 



1 sin [(A - X^)T/2] 



(A-A^) 



(42) 



This is the well-known sinc-function, which in the limit 
T — !> cxD is one of the definitions of the Dirac delta- 
function. The fact that the sinc-function Eq. ([^^ is 
sharply peaked around A = A^, and falls quickly when 
A — A-y > ^ gives rise to the notion of the resonance: 
substituted back into the integral over A in Eq. ([1T|) : 
the sinc-function suppresses all the values of the spec- 
tral density function a;Q,'y(A) outside this region. Thus, 
we observe that stretching the limits of integration over 
r' in Eq. (|41|1 to infinity implies the replacement of the 
actual sinc-function Eq. (|^^ by the Dirac delta-function 
i5(A — Xry). The error induced by this approximation is 



7 



legitimate to neglect for t = T/2 sufficiently large, and 
then Eq. (|4ip becomes the further approximated reduced 
master equation 



-27r 
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Transforming the last expression back to the Schrodinger 
representation we regain in matrix notation Eqs. p4l35p . 
except that instead of the constant density matrix lo we 
have the "resonant spectral density matrix" 



.Res 



N 

E 

a, 7 — 1 



(44) 



whose elements are given by cJq^(A) Eq. p6| eval- 
uated at the discrete eigenvalues A = {A^} of Hq, 
where the subscript 7 e [l,iV] corresponds to the 
column index. Thus, for large enough t = T/2 the 
contribution of the spectral density over the sharp 
resonances is very much enhanced, the structure of 
the rest of the continuous spectrum of Hq is quite 
unimportant for the reduced dynamics, and the error 
involved in the approximate Markovian Eq. (|43p is small. 

From the above considerations it is clear that the 
resonance-Markovian assumption fails if: (i) the interac- 
tion is strong, and (ii) uj{X) vanishes or undergoes signif- 
icant changes near the resonant values A ^ {A^}. Nev- 
ertheless, the Markovian approximation is usually very 
good and suits a large variety of natural environments. 
To see why, let us restore the Plank constant h, taken so 
far to be equal to unity. Doing so shows [ij] , that the re- 
gion of the non-negligible values of the sine-function, i.e. 
the resonance area, is equal to Since, for some rea- 
sonable T, the delta function approximation of the sinc- 
function is clearly a good one, the smoothness require- 
ment of Lu{X) refers to a finite number of small regions 
in the spectrum. The requirements of the approximation 
in Eq. (|43| are then not highly restrictive. Therefore, 
it would be experimentally difficult to detect any devi- 
ation from the semigroup evolution, which occurs only 
on "atomic" time-scales. In order to do this, one may 
need to deliberately destroy the smoothness of the envi- 
ronmental spectrum (by careful choice of the boundary 
conditions), which on such small scales might be difficult 
[23I nil- If ^(X) is changing rapidly (or vanishes) near 
the resonances, the Markovian approximation would be 
invalidated and the lifetime of the Zeno effect would be 
lengthened. 



gator U'^'^'^{t,tQ) expresses the fact that the dynamics of 
the open system may depend on its history. Yet, semi- 
group evolution for the open system can be achieved, 
if due to some special circumstances U^'^'^{t,to) may 
reduce to a dynamics approximated by the form of 
Eq. Since, in contrast to the closed system, the 

energy of the open system need not be conserved, the 
reduced evolution generator is not necessarily Hermitian 
and may have effective complex eigenvalues responsible 
for the decay. 

As mentioned, Eq. ([29]) may reduce to a Markovian 
form either exactly or approximately. The first case 
occurs when the spectral density matrix a;(A) Eq. ([26)) 
is constant and unbounded. The second, much more 
realistic, yields the desired effect approximately through 
a combination of the second order weak coupling per- 
turbation and resonances. In either case, the dramatic 
mutation of the time dependent spectral correlation 
matrix a{t) Eq. (|3T|) into a delta-correlated operator, 
leads to the equation of the form of Eq. (p4| . which 
independently of the dimension of the reduced subspace 
results in a semigroup solution of the form of Eq. psp . 
and thus, permits a semigroup decay also for the 
many channel problem, as shown in our analysis of the 
Wigner-Weisskopf method. 

To clarify the impact of the Markovian assumption 
on the pole approximation approach we first review 
its basics for the idealized Markovian limit. Note that 
in case the continuum spectrum is unbounded, there 
is no branch point, and the integration path C in the 
fundamental equation (|8]) goes above the whole real 
line. An immediate consequence of this fact is that the 



FIG. 3: The inverse Laplace transform contour C of Eq. ((8| in 
case of the unbounded perfect Markovian spectral coupling is 
modified into the contour Ci . The poles are denoted schemat- 
ically by the stars. 



WIGNER-WEISSKOPF POLE APPROXIMATION 
IN MARKOVIAN LIMIT 

In the preceding section we saw that the general 
integro-differential equation ()29l) for the reduced propa- 



modified integration path Ci of the background term 
J^^ R^^{z)e~'^^*dz in Eq. (jlOp may be taken parallel to 
the real axis Im 2 = 0, and if there is no obstacle to 
"dragging" it down to Im z — > —00, the background 
contribution vanishes for all positive times. Further, 
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we focus on the contributions of the pole terms. Recall 
Eq. ([7]) and note, that assuming this reduced equation 
is put in the Markovian form, implies that the evo- 
lution generator W{z) appearing in the denominator 
of Eq. ([7]) (similarly to the total system propagator 



Uiz) = Ct 



= i/{z — H)) is z-independent, even 



though not essentially Hcrmitian. The z-independence 
of W{z) automatically implies the z-independence of its 
eigenvectors and eigenvalues defined in Eq. (|13p . and 
hence, guarantees the orthogonality of the projectors 
Qaj{zj), Qa^izk) if 7^ uja^ for any Zj,Zk. Therefore, 
the idealized Markovian limit allows pure semigroup 
evolution with no Zcno effect [2^ independently of the 
position or the number of the poles. 

To prove our last conclusion explicitly we establish the 
connection between the reduced generator W^^{z) and 
the spectral correlation matrix in the lower half of the 
Laplace plane. Combining Eq. ()9IA.7IA.8|) we find 



ia{z). 



(45) 



Generally, the spectral correlation matrix a{t) (Eq. pip ), 
standing for the memory kernel of the integro-differential 
equation ([29]), may be any spread function of time. The 
essence of the Markovian limit consists of collapsing (ei- 
ther exactly or approximately) this memory-kernel into a 
delta-correlated matrix. For the "artificially" unbounded 
spectrum A G (—00, 00) we have 



a 



and hence, 



Markov / 



(t) / dAw(A)e 



-iXt 



2TrojS{t), (46) 



■^Mar/coo (^-j = i2TrLJ / dtS{t)e''' - = -TTiuJ. 

J —00 

(47) 

where the spectral density matrix cj is a constant. Sub- 
stituting this back into Ec^. yields 



' (z) = H^'"'^ + ia^^'^^'^"^ (z) = const. (48) 



What is left is to return to the original formulation 
of the Wigner-Wcisskopf theory with the semi-bounded 
continuous spectrum and explain how the effect of weak 
coupling and resonance explain the high precision of 
semigroup evolution either for the single or for the 
many channel decay. In reaching the above result, it 
was assumed w(A) is supported in (— cxd,oo). We prove 
in Appendix A that h^^{z) may have zero determinant 
for some values of z in the lower half plane. The 
leading behavior of the reduced evolution is therefore 
dictated only by a finite number of the small regions 
corresponding to the poles of R^^{z), where the contri- 
bution of the spectral density matrix w(A) is strongly 
enhanced, in agreement with resonance-Markovian 
approximation. In such a case, the integral over A 



in Eq. (|A.9[) contributes primarily on the set of A's 
close to the real part of the poles, whereas a fictitious 
extension of the integration to —00 will not change 
appreciably the value of the integral. Hence, it is 
straightforward to repeat the derivations leading to 
the conclusion on the simultaneous vanishing of the 
background contribution and the orthogonality of the 
pole term projectors. The only difference is that instead 
of the constant spectral density matrix ui in Eqs. ()46I47|) 
we shall have the resonant A-independent matrix cj^*^* 
given by Eq. To conclude, the result Eq. (|i5| is 

quite general in case the resonances are fairly sharp (i.e., 
the poles in the second sheet are close to the real line). 
Even though not strictly exact, this Markovian limit 
approximation should be valid for many experimental 
conditions, and becomes even better for weaker coupling. 



SUMMARY 

By associating the spectral weights of the Wigner- 
Weisskopf model for unstable system decay with the sta- 
tistical properties of coupling to an environment, we were 
able to characterize the reduced evolution in the sub- 
space of unstable states in terms of the spectral correla- 
tion matrix. Exploiting the properties of the latter, we 
showed that the Markovian limit distribution is sufficient 
to account for semigroup behavior for an arbitrary num- 
ber of the decay channels, observed in experiment, such 
as in 0. Besides the important impact of the spectral 
correlation matrix notion for understanding the reduced 
evolution of the multi-channel decays, such as the two- 
channel K- meson decay 18 , 3] , it is clear that the asso- 
ciation of the coupling to a spectral continuum with an 
environment may be similarly useful for any analogously 
modeled theory. 



ACKNOWLEDGEMENTS 

We thank Dr. Y. Strauss for discussion of analogous 
results obtained by Nagy-Foias-Kolmogorov dilation [25| , 
in which the quantum mechanical problem is imbedded 
in a larger space as in the work of Maassen [i^, who 
constructs a Fock space resulting in the generation of a 
"noisy environment" of the type we have discussed. 



APPENDIX A. SECOND SHEET PROPERTIES 
OF R(z) 

In this appendix we discuss the structure of the re- 
duced propagator in the complex Laplace plane. We 
specifically show that the propagator has poles only in 
the lower half plane near the real axis and also provide 
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the background needed for establishing the connection 
between the reduced system evolution generator and the 
spectral density matrix Eq. (|^5|) . Following the well- 
known procedure let us rewrite Eq. ([7]) as 



R{z) 
Here the operator 

h{z) ^z- H^''^ 



h{zy 



dA 



■z-A' 



(A.l) 



(A.2) 



with Hq"'- and given by Eqs. (|30I32I) respectively, 

is found straightforwardly from the projection of the to- 
tal system propagator on the discrete subspace. The last 
term of Eq. (jA.2p can be recognized as the Laplace trans- 
form (for Imz > 0) of the spectral correlation matrix a{t) 
Eq. dSI]): 



ia{z) 



dX 



c.(A) 



i / dta(t)e'' 

'0 



dt / dAa;(A)e^ 



-iXt ^izt 



(A.3) 



To obtain a singularity in R{z) Eq. (jA.l[) we need that 
the determinant of h{z) Eq. (jA.2[) vanishes. These non- 
Hermitian matrices can be put into Jordan canonical 
form with a unitary transformation, with eigenvalues 
along the diagonal. A vanishing eigenvalue implies a van- 
ishing determinant. We therefore can ask whether h{z) 
has a vanishing eigenvalue, i.e. 



h{z)\x{z))K= z - Hl^,'" - I |xW)r = 0, 

. (^-^^ 

where |x(-2))r is the right eigenvector of h{z) with the 
presumed zero eigenvalue. Taking the scalar product 
with eigenvector from the left yields 

z\x{z))? -i.{x{z)\HI^''\x{z))k~ 



A 



i.e., Im z times a positive quantity. This can never be 
zero for Im z on the first Riemann sheet. We are 
therefore required to go to the second sheet: 



R'\z) = 



1 



= z~ H'o"'^ - ia{z). 



(A.7) 



(A., 



Here R^^{z) and h{zY^ are the smooth continuations of 
R{z) and h{z) into the lower half plane and 



ia{z) = i dt dXuj{X)e 



-iXtAzt 



2muj{z). (A.9) 



is the analytic continuation of the spectral correlation 
matrix a{z) Eq. (|A.3p . where lo{z) is the analytic contin- 
uation of w(A) into the lower half plane. The condition 
for a vanishing eigenvalue Eq. (|A.5[) now reads 

^^ L(x(^)|^(A)|x(z))J^ , o ■ / / M ^ M ^ ^\ n 
dA +2Tni,{x{z)\u{z)\x{z)) = 0, 

Z — A 

Taking the imaginary part as before gives 

L(x(^)k(A)|xW)i^ 



(A.IO) 



Imz<^ \x{z))V + / dA 



|z-A|^ 



2'Kii^{x{z)\uj{z)\x{z)) I 



0, 



(A.ll) 



so that there may be a solution for Imz negative, since 
a; (A) is positive definite. It is easy to see that the 
expectation of a;(A) is the sum of absolute squares; the 
analytic continuation of lo{X) to the lower half plane for 
small imaginary part of z, enough to reach a resonance 
pole, is assumed to remain approximately real since it is 
smoothly connected to real values of a; (A), for A on the 
real line. In the Markovian limit this function is taken 
to approach a constant. 



- / dAi:MfMAMf))£ ^ ^ (A.5) 
J z- X 

In general, for the convergence of the integral in the last 
expression, we must assume that ui{X) decreases better 
than an, no matter how small, inverse power of A since 
the integral contains ^ for very large A. Now we 
consider the imaginary part of Eq. (jA.Sp . Since H^'^'^ is 
Hermitian, its expectation value in the state Ixiz))^ is 
real, and we are left with 

Imz{|x(.))P+/dAiiMf)M^W^ (A.6) 
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